Whole genome sequencing based analysis of inflammation biomarkers in the Trans-Omics for Precision Medicine (TOPMed) consortium.

Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38 465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program (with varying sample size by trait, where the minimum sample size was n = 737 for MMP-1). We identified 22 distinct single-variant associations across 6 traits-E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin-that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

Whole Genome Sequencing Based Analysis of Inflammation Biomarkers in the Trans-Omics for Precision Medicine (TOPMed) Consortium.

Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38,465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program. We identified 22 distinct single-variant associations across 6 traits - E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin - that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

Characterization of POT1 tumor predisposition syndrome: Tumor prevalence in a clinically diverse hereditary cancer cohort.

PURPOSE: Germline variants in POT1 have been implicated in predisposition to melanoma, sarcoma, and glioma in limited studies. Here, we determine the prevalence of cancer types in individuals with POT1 pathogenic variants (PVs) undergoing multigene panel testing (MGPT) for a broad variety of cancer indications. METHODS: We performed a retrospective review of data provided on clinical documents from individuals with POT1 PVs identified via MGPT over a 5-year period. Tumor prevalence in POT1 PV heterozygotes was compared with MGPT-negative wild-type (WT) controls using χ2 test. RESULTS: POT1 PVs were identified in 227 individuals. POT1 PV and WT (n = 13,315) cohorts had a similar proportion of reported tumors (69.6% and 69.2%, respectively); however, POT1 PV heterozygotes were more likely to be diagnosed with multiple tumors (18.9% vs 8.7%; P < .001). Compared with POT1 WT, we identified a significant increase in melanoma (odds ratio 7.03; 95% CI 4.7-10.5; P < .001) and sarcoma (odds ratio 6.6; 95% CI 3.1-13.9; P < .001). CONCLUSION: This analysis of the largest POT1 PV cohort to date validates the inclusion of POT1 in hereditary cancer MGPT and has the potential to impact clinical management recommendations, particularly for patients and families at risk for melanoma and sarcoma.

Multi-ethnic GWAS and fine-mapping of glycaemic traits identify novel loci in the PAGE Study.

AIMS/HYPOTHESIS: Type 2 diabetes is a growing global public health challenge. Investigating quantitative traits, including fasting glucose, fasting insulin and HbA1c, that serve as early markers of type 2 diabetes progression may lead to a deeper understanding of the genetic aetiology of type 2 diabetes development. Previous genome-wide association studies (GWAS) have identified over 500 loci associated with type 2 diabetes, glycaemic traits and insulin-related traits. However, most of these findings were based only on populations of European ancestry. To address this research gap, we examined the genetic basis of fasting glucose, fasting insulin and HbA1c in participants of the diverse Population Architecture using Genomics and Epidemiology (PAGE) Study. METHODS: We conducted a GWAS of fasting glucose (n = 52,267), fasting insulin (n = 48,395) and HbA1c (n = 23,357) in participants without diabetes from the diverse PAGE Study (23% self-reported African American, 46% Hispanic/Latino, 40% European, 4% Asian, 3% Native Hawaiian, 0.8% Native American), performing transethnic and population-specific GWAS meta-analyses, followed by fine-mapping to identify and characterise novel loci and independent secondary signals in known loci. RESULTS: Four novel associations were identified (p < 5 × 10-9), including three loci associated with fasting insulin, and a novel, low-frequency African American-specific locus associated with fasting glucose. Additionally, seven secondary signals were identified, including novel independent secondary signals for fasting glucose at the known GCK locus and for fasting insulin at the known PPP1R3B locus in transethnic meta-analysis. CONCLUSIONS/INTERPRETATION: Our findings provide new insights into the genetic architecture of glycaemic traits and highlight the continued importance of conducting genetic studies in diverse populations. DATA AVAILABILITY: Full summary statistics from each of the population-specific and transethnic results are available at NHGRI-EBI GWAS catalog ( https://www.ebi.ac.uk/gwas/downloads/summary-statistics ).

Genetic discovery and risk characterization in type 2 diabetes across diverse populations.

Genomic discovery and characterization of risk loci for type 2 diabetes (T2D) have been conducted primarily in individuals of European ancestry. We conducted a multiethnic genome-wide association study of T2D among 53,102 cases and 193,679 control subjects from African, Hispanic, Asian, Native Hawaiian, and European population groups in the Population Architecture Genomics and Epidemiology (PAGE) and Diabetes Genetics Replication and Meta-analysis (DIAGRAM) Consortia. In individuals of African ancestry, we discovered a risk variant in the TGFB1 gene (rs11466334, risk allele frequency (RAF) = 6.8%, odds ratio [OR] = 1.27, p = 2.06 × 10-8), which replicated in independent studies of African ancestry (p = 6.26 × 10-23). We identified a multiethnic risk variant in the BACE2 gene (rs13052926, RAF = 14.1%, OR = 1.08, p = 5.75 × 10-9), which also replicated in independent studies (p = 3.45 × 10-4). We also observed a significant difference in the performance of a multiethnic genetic risk score (GRS) across population groups (pheterogeneity = 3.85 × 10-20). Comparing individuals in the top GRS risk category (40%-60%), the OR was highest in Asians (OR = 3.08) and European (OR = 2.94) ancestry populations, followed by Hispanic (OR = 2.39), Native Hawaiian (OR = 2.02), and African ancestry (OR = 1.57) populations. These findings underscore the importance of genetic discovery and risk characterization in diverse populations and the urgent need to further increase representation of non-European ancestry individuals in genetics research to improve genetic-based risk prediction across populations.

Whole-genome sequencing in diverse subjects identifies genetic correlates of leukocyte traits: The NHLBI TOPMed program.

Many common and rare variants associated with hematologic traits have been discovered through imputation on large-scale reference panels. However, the majority of genome-wide association studies (GWASs) have been conducted in Europeans, and determining causal variants has proved challenging. We performed a GWAS of total leukocyte, neutrophil, lymphocyte, monocyte, eosinophil, and basophil counts generated from 109,563,748 variants in the autosomes and the X chromosome in the Trans-Omics for Precision Medicine (TOPMed) program, which included data from 61,802 individuals of diverse ancestry. We discovered and replicated 7 leukocyte trait associations, including (1) the association between a chromosome X, pseudo-autosomal region (PAR), noncoding variant located between cytokine receptor genes (CSF2RA and CLRF2) and lower eosinophil count; and (2) associations between single variants found predominantly among African Americans at the S1PR3 (9q22.1) and HBB (11p15.4) loci and monocyte and lymphocyte counts, respectively. We further provide evidence indicating that the newly discovered eosinophil-lowering chromosome X PAR variant might be associated with reduced susceptibility to common allergic diseases such as atopic dermatitis and asthma. Additionally, we found a burden of very rare FLT3 (13q12.2) variants associated with monocyte counts. Together, these results emphasize the utility of whole-genome sequencing in diverse samples in identifying associations missed by European-ancestry-driven GWASs.

Genome-wide association study of pancreatic fat: The Multiethnic Cohort Adiposity Phenotype Study.

Several studies have found associations between higher pancreatic fat content and adverse health outcomes, such as diabetes and the metabolic syndrome, but investigations into the genetic contributions to pancreatic fat are limited. This genome-wide association study, comprised of 804 participants with MRI-assessed pancreatic fat measurements, was conducted in the ethnically diverse Multiethnic Cohort-Adiposity Phenotype Study (MEC-APS). Two genetic variants reaching genome-wide significance, rs73449607 on chromosome 13q21.2 (Beta = -0.67, P = 4.50x10-8) and rs7996760 on chromosome 6q14 (Beta = -0.90, P = 4.91x10-8) were associated with percent pancreatic fat on the log scale. Rs73449607 was most common in the African American population (13%) and rs79967607 was most common in the European American population (6%). Rs73449607 was also associated with lower risk of type 2 diabetes (OR = 0.95, 95% CI = 0.89-1.00, P = 0.047) in the Population Architecture Genomics and Epidemiology (PAGE) Study and the DIAbetes Genetics Replication and Meta-analysis (DIAGRAM), which included substantial numbers of non-European ancestry participants (53,102 cases and 193,679 controls). Rs73449607 is located in an intergenic region between GSX1 and PLUTO, and rs79967607 is in intron 1 of EPM2A. PLUTO, a lncRNA, regulates transcription of an adjacent gene, PDX1, that controls beta-cell function in the mature pancreas, and EPM2A encodes the protein laforin, which plays a critical role in regulating glycogen production. If validated, these variants may suggest a genetic component for pancreatic fat and a common etiologic link between pancreatic fat and type 2 diabetes.

Risk of breast cancer and prediagnostic urinary excretion of bisphenol A, triclosan and parabens: The Multiethnic Cohort Study.

Exposure to bisphenol A (BPA), triclosan and parabens is widespread but their impact on breast cancer risk remains unclear. This nested case-control study investigated endocrine-disrupting chemicals (EDCs) and breast cancer risk within the Multiethnic Cohort (MEC). We measured prediagnostic urinary BPA, triclosan and parabens in 1032 postmenopausal women with breast cancer (48 African American, 77 Latino, 155 Native Hawaiian, 478 Japanese American and 274 White) and 1030 individually matched controls, using a sensitive and validated liquid chromatography mass spectrometry assay. Conditional logistic regression was used to examine risk with these EDCs with adjustment for creatinine and potential confounders. In all women, breast cancer risk was not associated with BPA (Ptrend  = 0.53) and was inversely associated with triclosan (ORT3 vs T1  = 0.83, 95% CI: 0.66-1.04, Ptrend  = 0.045) and total parabens (ORT3 vs T1  = 0.77, 95% CI: 0.62-0.97, Ptrend  = 0.03). While risk of hormone receptor positive (HR+) cancer was 20% to 23% lower among women in the upper two tertiles of paraben exposure (Ptrend  = 0.02), risk of HR negative (HR-) was reduced 27% but only among those in the upper tertile of exposure. Although risk associations did not differ significantly by ethnicity or by body mass index (BMI), the inverse association with triclosan was observed mainly among overweight/obese women (ORT3 vs T1  = 0.76, 95% CI: 0.56-1.02, Ptrend  = 0.02). In summary, breast cancer risk in a multiethnic population was unrelated to BPA and was weakly inversely associated with triclosan and paraben exposures. Studies with multiple urine samples collected before breast cancer diagnosis are needed to further investigate these EDCs and breast cancer risk.

Urinary phthalate exposures and risk of breast cancer: the Multiethnic Cohort study.

BACKGROUND: The epidemiologic evidence from observational studies on breast cancer risk and phthalates, endocrine disrupting chemicals, has been inconsistent. In the only previous study based on pre-diagnostic urinary phthalates and risk of breast cancer, results were null in mostly white women. METHODS: We examined the association between pre-diagnostic urinary phthalates and breast cancer in a nested case-control study within the Multiethnic Cohort (MEC) study, presenting the first data from five major racial/ethnic groups in the USA. We measured 10 phthalate metabolites and phthalic acid, using a sensitive liquid chromatography mass spectrometry assay on 1032 women with breast cancer (48 African Americans, 77 Latinos, 155 Native Hawaiians, 478 Japanese Americans, and 274 Whites) and 1030 matched controls. Conditional logistic regression was used to examine risk with individual metabolites and ratios of primary (MEHP, mono-2-ethylhexyl-phthalate) to secondary (MEHHP, mono(2-ethyl-5-hydroxyhexyl); MEOHP, mono(2-ethyl-5-oxohexy)) metabolites of di-2-ethylhexyl phthalate (DEHP), a widely used plasticizer. In addition, we investigated risk associations with high (∑HMWP) and low molecular weight (∑LMWP) phthalates, as well as total phthalates which included high and low molecular weight phthalates with phthalic acid (∑LMHMPA) or without phthalic acid in molar ratios (∑LMHMmolar) and adjusted for creatinine and potential confounders. RESULTS: Among all women, breast cancer risk was higher for those in tertile 2 and tertile 3 of primary to secondary metabolites of DEHP (MEHP/(MEHHP + MEOHP)) in comparison to those in tertile 1; the respective odds ratios were 1.32 (95% CI 1.04-1.68) and 1.26 (95% CI 0.96-1.66) (Ptrend = 0.05). Risk among Native Hawaiian women increased with exposures to eight of ten individual phthalates and total phthalates (∑LMHMPA ORT3 vs T1 = 2.66, 95% CI 1.39-5.12, Ptrend = 0.001). In analysis by hormone receptor (HR) status, exposure above the median of ∑LMWP was associated with an increased risk of HR-positive breast cancer (OR = 1.30, 95% CI 1.05-1.60) while above the median exposure to phthalic acid was associated with an increased risk of HR-negative breast cancer (ORabove vs below median = 1.59, 95% CI 1.01-2.48). CONCLUSIONS: Further investigations of suggestive associations of elevated breast cancer risk with higher ratios of primary to secondary metabolites of DEHP, and differences in risk patterns by race/ethnicity and HR status are warranted.

Allelic Heterogeneity at the CRP Locus Identified by Whole-Genome Sequencing in Multi-ancestry Cohorts.

Whole-genome sequencing (WGS) can improve assessment of low-frequency and rare variants, particularly in non-European populations that have been underrepresented in existing genomic studies. The genetic determinants of C-reactive protein (CRP), a biomarker of chronic inflammation, have been extensively studied, with existing genome-wide association studies (GWASs) conducted in >200,000 individuals of European ancestry. In order to discover novel loci associated with CRP levels, we examined a multi-ancestry population (n = 23,279) with WGS (∼38× coverage) from the Trans-Omics for Precision Medicine (TOPMed) program. We found evidence for eight distinct associations at the CRP locus, including two variants that have not been identified previously (rs11265259 and rs181704186), both of which are non-coding and more common in individuals of African ancestry (∼10% and ∼1% minor allele frequency, respectively, and rare or monomorphic in 1000 Genomes populations of East Asian, South Asian, and European ancestry). We show that the minor (G) allele of rs181704186 is associated with lower CRP levels and decreased transcriptional activity and protein binding in vitro, providing a plausible molecular mechanism for this African ancestry-specific signal. The individuals homozygous for rs181704186-G have a mean CRP level of 0.23 mg/L, in contrast to individuals heterozygous for rs181704186 with mean CRP of 2.97 mg/L and major allele homozygotes with mean CRP of 4.11 mg/L. This study demonstrates the utility of WGS in multi-ethnic populations to drive discovery of complex trait associations of large effect and to identify functional alleles in noncoding regulatory regions.

Whole genome sequence association with E-selectin levels reveals loss-of-function variant in African Americans.

E-selectin mediates the rolling of circulating leukocytes during inflammatory processes. Previous genome-wide association studies in European and Asian individuals have identified the ABO locus associated with E-selectin levels. Using Trans-Omics for Precision Medicine whole genome sequencing data in 2249 African Americans (AAs) from the Jackson Heart Study, we examined genome-wide associations with soluble E-selectin levels. In addition to replicating known signals at ABO, we identified a novel association of a common loss-of-function, missense variant in Fucosyltransferase 6 (FUT6; rs17855739,p.Glu274Lys, P = 9.02 × 10-24) with higher soluble E-selectin levels. This variant is considerably more common in populations of African ancestry compared to non-African ancestry populations. We replicated the association of FUT6 p.Glu274Lys with higher soluble E-selectin in an independent population of 748 AAs from the Women's Health Initiative and identified an additional pleiotropic association with vitamin B12 levels. Despite the broad role of both selectins and fucosyltransferases in various inflammatory, immune and cancer-related processes, we were unable to identify any additional disease associations of the FUT6 p.Glu274Lys variant in an electronic medical record-based phenome-wide association scan of over 9000 AAs.

Lipoprotein-associated phospholipase A2 and risk of incident peripheral arterial disease: Findings from The Atherosclerosis Risk in Communities study (ARIC).

BACKGROUND AND AIMS: Results from prospective studies evaluating the relationship between elevated lipoprotein-associated phospholipase A2 (Lp-PLA2) activity and incident peripheral arterial disease (PAD) have been mixed. We investigated whether higher Lp-PLA2 levels are associated with increased risk of incident PAD and whether PLA2G7 gene variants, which result in lower Lp-PLA2 levels, are associated with reduced risk of incident PAD. METHODS: Our analysis included 9922 participants (56% female; 21% African-American; mean age 63 years) without baseline PAD at ARIC Visit 4 (1996-1998), who had Lp-PLA2 activity measured and were subsequently followed for the development of PAD, defined by occurrence of a PAD-related hospitalization, through 2012. Cox proportional hazard models were performed to determine the association of Lp-PLA2 levels and PLA2G7 gene variants with incident PAD. RESULTS: During a median follow-up of 14.9 years, we identified 756 incident cases of PAD. In analyses adjusting for age, race, and sex, each standard deviation increment in Lp-PLA2 activity (62 nmol/ml/min) was associated with a higher risk of developing PAD (hazard ratio (HR) 1.17; 95% confidence interval (CI) 1.09, 1.26). This association remained significant after additional adjustment for risk factors, other cardiovascular disease, and medication use, but was strongly attenuated (HR: 1.09; 95% CI 1.00, 1.20). PLA2G7 variants were not associated with a lower risk of PAD in both white carriers (HR: 1.21; 95% CI: 0.17-8.56) and African-American carriers (HR: 0.83; 95% CI: 0.41-1.67), although statistical power was quite limited for this analysis, particularly in whites. CONCLUSIONS: While higher Lp-PLA2 activity was associated with an increased risk for incident PAD, it is likely a risk marker largely represented by traditional risk factors.

Genetic invalidation of Lp-PLA2 as a therapeutic target: Large-scale study of five functional Lp-PLA2-lowering alleles.

Aims Darapladib, a potent inhibitor of lipoprotein-associated phospholipase A2 (Lp-PLA2), has not reduced risk of cardiovascular disease outcomes in recent randomized trials. We aimed to test whether Lp-PLA2 enzyme activity is causally relevant to coronary heart disease. Methods In 72,657 patients with coronary heart disease and 110,218 controls in 23 epidemiological studies, we genotyped five functional variants: four rare loss-of-function mutations (c.109+2T > C (rs142974898), Arg82His (rs144983904), Val279Phe (rs76863441), Gln287Ter (rs140020965)) and one common modest-impact variant (Val379Ala (rs1051931)) in PLA2G7, the gene encoding Lp-PLA2. We supplemented de-novo genotyping with information on a further 45,823 coronary heart disease patients and 88,680 controls in publicly available databases and other previous studies. We conducted a systematic review of randomized trials to compare effects of darapladib treatment on soluble Lp-PLA2 activity, conventional cardiovascular risk factors, and coronary heart disease risk with corresponding effects of Lp-PLA2-lowering alleles. Results Lp-PLA2 activity was decreased by 64% ( p = 2.4 × 10-25) with carriage of any of the four loss-of-function variants, by 45% ( p < 10-300) for every allele inherited at Val279Phe, and by 2.7% ( p = 1.9 × 10-12) for every allele inherited at Val379Ala. Darapladib 160 mg once-daily reduced Lp-PLA2 activity by 65% ( p < 10-300). Causal risk ratios for coronary heart disease per 65% lower Lp-PLA2 activity were: 0.95 (0.88-1.03) with Val279Phe; 0.92 (0.74-1.16) with carriage of any loss-of-function variant; 1.01 (0.68-1.51) with Val379Ala; and 0.95 (0.89-1.02) with darapladib treatment. Conclusions In a large-scale human genetic study, none of a series of Lp-PLA2-lowering alleles was related to coronary heart disease risk, suggesting that Lp-PLA2 is unlikely to be a causal risk factor.

Whole-exome sequencing reveals an inherited R566X mutation of the epithelial sodium channel ß-subunit in a case of early-onset phenotype of Liddle syndrome.

To comprehensively evaluate a European-American child with severe hypertension, whole-exome sequencing (WES) was performed on the child and parents, which identified causal variation of the proband's early-onset disease. The proband's hypertension was resistant to treatment, requiring a multiple drug regimen including amiloride, spironolactone, and hydrochlorothiazide. We suspected a monogenic form of hypertension because of the persistent hypokalemia with low plasma levels of renin and aldosterone. To address this, we focused on rare functional variants and indels, and performed gene-based tests incorporating linkage scores and allele frequency and filtered on deleterious functional mutations. Drawing upon clinical presentation, 27 genes were selected evidenced to cause monogenic hypertension and matched to the gene-based results. This resulted in the identification of a stop-gain mutation in an epithelial sodium channel (ENaC), SCNN1B, an established Liddle syndrome gene, shared by the child and her father. Interestingly, the father also harbored a missense mutation (p.Trp552Arg) in the α-subunit of the ENaC trimer, SCNN1A, possibly pointing to pseudohypoaldosteronism type I. This case is unique in that we present the early-onset disease and treatment response caused by a canonical stop-gain mutation (p.Arg566*) as well as ENaC digenic hits in the father, emphasizing the utility of WES informing precision medicine.

Whole-Exome Sequencing Identifies Loci Associated with Blood Cell Traits and Reveals a Role for Alternative GFI1B Splice Variants in Human Hematopoiesis.

Circulating blood cell counts and indices are important indicators of hematopoietic function and a number of clinical parameters, such as blood oxygen-carrying capacity, inflammation, and hemostasis. By performing whole-exome sequence association analyses of hematologic quantitative traits in 15,459 community-dwelling individuals, followed by in silico replication in up to 52,024 independent samples, we identified two previously undescribed coding variants associated with lower platelet count: a common missense variant in CPS1 (rs1047891, MAF = 0.33, discovery + replication p = 6.38 × 10(-10)) and a rare synonymous variant in GFI1B (rs150813342, MAF = 0.009, discovery + replication p = 1.79 × 10(-27)). By performing CRISPR/Cas9 genome editing in hematopoietic cell lines and follow-up targeted knockdown experiments in primary human hematopoietic stem and progenitor cells, we demonstrate an alternative splicing mechanism by which the GFI1B rs150813342 variant suppresses formation of a GFI1B isoform that preferentially promotes megakaryocyte differentiation and platelet production. These results demonstrate how unbiased studies of natural variation in blood cell traits can provide insight into the regulation of human hematopoiesis.

Association of Exome Sequences With Cardiovascular Traits Among Blacks in the Jackson Heart Study.

BACKGROUND: The correlation of null alleles with human phenotypes can provide insight into gene function in humans. In individuals of African ancestry, we set out to identify null and damaging missense variants, and test these variants for association with a range of cardiovascular phenotypes. METHODS AND RESULTS: We performed whole-exome sequencing in 3223 black individuals from the Jackson Heart Study and found a total of 729 666 variant sites with minor allele frequency <5%, including 17 263 null variants and 49 929 missense variants predicted to be damaging by in silico algorithms. We tested null and damaging missense variants within each gene for association with 36 cardiovascular traits. We found 3 associations that met our prespecified level of significance (α=1.1×10(-7)). Null and damaging missense variants in PCSK9 were associated with 36 mg/dL lower low-density lipoprotein cholesterol (P=3×10(-21)). Three individuals in their 50s with complete PCSK9 deficiency (each compound heterozygote for PCSK9 p.Y142X and p.C679X) were identified, with one having a coronary artery calcification score in the 83rd percentile despite a low-density lipoprotein cholesterol of 32 mg/dL. A damaging missense variant in HBQ1 (p.G52A) was associated with a 2 pg/cell lower mean corpuscular hemoglobin (P=9×10(-13)) and rare damaging missense variants in VPS13A with higher red blood cell distribution width (P=9.9×10(-8)). CONCLUSIONS: A limited number of null/damaging alleles with a large effect on cardiovascular traits were detectable in ≈3000 black individuals.

Platelet-Related Variants Identified by Exomechip Meta-analysis in 157,293 Individuals.

Platelet production, maintenance, and clearance are tightly controlled processes indicative of platelets' important roles in hemostasis and thrombosis. Platelets are common targets for primary and secondary prevention of several conditions. They are monitored clinically by complete blood counts, specifically with measurements of platelet count (PLT) and mean platelet volume (MPV). Identifying genetic effects on PLT and MPV can provide mechanistic insights into platelet biology and their role in disease. Therefore, we formed the Blood Cell Consortium (BCX) to perform a large-scale meta-analysis of Exomechip association results for PLT and MPV in 157,293 and 57,617 individuals, respectively. Using the low-frequency/rare coding variant-enriched Exomechip genotyping array, we sought to identify genetic variants associated with PLT and MPV. In addition to confirming 47 known PLT and 20 known MPV associations, we identified 32 PLT and 18 MPV associations not previously observed in the literature across the allele frequency spectrum, including rare large effect (FCER1A), low-frequency (IQGAP2, MAP1A, LY75), and common (ZMIZ2, SMG6, PEAR1, ARFGAP3/PACSIN2) variants. Several variants associated with PLT/MPV (PEAR1, MRVI1, PTGES3) were also associated with platelet reactivity. In concurrent BCX analyses, there was overlap of platelet-associated variants with red (MAP1A, TMPRSS6, ZMIZ2) and white (PEAR1, ZMIZ2, LY75) blood cell traits, suggesting common regulatory pathways with shared genetic architecture among these hematopoietic lineages. Our large-scale Exomechip analyses identified previously undocumented associations with platelet traits and further indicate that several complex quantitative hematological, lipid, and cardiovascular traits share genetic factors.

Large-Scale Exome-wide Association Analysis Identifies Loci for White Blood Cell Traits and Pleiotropy with Immune-Mediated Diseases.

White blood cells play diverse roles in innate and adaptive immunity. Genetic association analyses of phenotypic variation in circulating white blood cell (WBC) counts from large samples of otherwise healthy individuals can provide insights into genes and biologic pathways involved in production, differentiation, or clearance of particular WBC lineages (myeloid, lymphoid) and also potentially inform the genetic basis of autoimmune, allergic, and blood diseases. We performed an exome array-based meta-analysis of total WBC and subtype counts (neutrophils, monocytes, lymphocytes, basophils, and eosinophils) in a multi-ancestry discovery and replication sample of âˆ¼157,622 individuals from 25 studies. We identified 16 common variants (8 of which were coding variants) associated with one or more WBC traits, the majority of which are pleiotropically associated with autoimmune diseases. Based on functional annotation, these loci included genes encoding surface markers of myeloid, lymphoid, or hematopoietic stem cell differentiation (CD69, CD33, CD87), transcription factors regulating lineage specification during hematopoiesis (ASXL1, IRF8, IKZF1, JMJD1C, ETS2-PSMG1), and molecules involved in neutrophil clearance/apoptosis (C10orf54, LTA), adhesion (TNXB), or centrosome and microtubule structure/function (KIF9, TUBD1). Together with recent reports of somatic ASXL1 mutations among individuals with idiopathic cytopenias or clonal hematopoiesis of undetermined significance, the identification of a common regulatory 3' UTR variant of ASXL1 suggests that both germline and somatic ASXL1 mutations contribute to lower blood counts in otherwise asymptomatic individuals. These association results shed light on genetic mechanisms that regulate circulating WBC counts and suggest a prominent shared genetic architecture with inflammatory and autoimmune diseases.